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ABSTRACT 

Context. There are several wide field galaxy and cluster surveys planned for the nearest future, e.g. BOSS, WFMOS, ADEPT, Hetdex, 
SPT, eROSITA. In the simplest approach one would analyze these independently, thus neglecting the extra information provided by 
the cluster-galaxy cross-pairs. 

Aims. In this paper we have focused on the possible synergy between these surveys by investigating the amount of information en- 
coded in the cross-pairs. 

Methods. We present a model for the cluster-galaxy cross-spectrum within the Halo Model framework. To assess the gain in perfor- 
mance due to inclusion of the cluster-galaxy cross-pairs we carry out a Fisher matrix analysis for a BOSS-like galaxy redshift survey 
targeting luminous red galaxies and a hypothetical mass-limited cluster redshift survey with a lower mass threshold of 1.7x lO'** /?"'Mq 
over the same volume. 

Results. On small scales cluster-galaxy cross-spectrum probes directly density profile of the halos, instead of the density profile con- 
volved with itself, as is the case for the galaxy power spectrum. Due to this different behavior, adding information from the cross-pairs 
helps to tighten constraints on the halo occupation distribution (e.g. a factor of ~ 2 compression of the error ellipses on the m'°'"-a 
plane) and offers an alternative mechanism compared with techniques that directly fit halo density profiles. By inclusion of the cross- 
pairs a factor of ~ 2 stronger constraints are obtained for erg, while the improvement for the dark energy figure-of-merit is somewhat 
weaker: an increase by a factor of 1.4. 

We have also written down the formalism for the case when only photometric redshifts are available for both the clusters and the 
galaxies. For the analysis of the photometric surveys the inclusion of the cluster-galaxy cross-pairs might be very beneficial since the 
photo-z errors for the clusters are usually significantly smaller than that for the typical galaxies. 

Key words. Cosmology: theory - large-scale structure of Universe - Galaxies: clusters: general ~ Galaxies: statistics - cosmological 
parameters 



1. Introduction 

With the advent of the modern wide field spectroscopic sky sur- 
veys, which have measured redshifts for hundreds of thousands 
of galaxies, e.g. 2dFGRSEl SDSSB we have gained a good pic- 
ture of the cosmic large-scale structure out to redshifts of ~ 0.5. 
From the measurements of the cosmi c microwave backgro und 
angular temperature fluctuations (e.g. iHinshaw et aP 1200 8*) we 
have obtained a good understanding of the nature and statistics 
of the small density fluctuations that through gravitational in- 
stability have evolved to the filament-void network of the mat- 
ter surrounding us at low redshifts. These measurements have 
proven that initial fluctuations are compatible of being adiabatic 
and Gaussian with roughly the scale-free initial spectrum. As 
the initial fluctuations can be described as a Gaussian random 
field the full information is contained in the two-point corre- 
lator of the field. Under the linear evolution of the initial fluc- 
tuation field, which is valid on large scales (k < 0.1 /iMpc"'), 
the Gaussian field remains Gaussian in nature, and as such ex- 
plains the central interest in measuring the clustering power 
spectrum (or its real space analog - two-point correlation func- 



tion) using the high quality data from the above-mentioned 
spectroscopic surveys. Data from the 2dFGRS and from the 
SDSS survey have both lead to the high precis ion measurement 
of t he low-redshift g a laxy power spect ra (Percival et al.ll200ll 
Teg market all llOOl ICole et al.1 llOOl [jegmarket al. 200^. 
Moreover, these datasets have proven to be powerful enough 
to enable one to detect theoretically expected small fluctua- 
tions (~ 5% relative amplitude) - baryonic acou stic oscilla- 
tions (B AO) - in the clusteri ng pow er spectra ( Eisenstein et alj 
2005; C ole et al.l |2p05; Hut sill2006al: iPadmana bhan et al.ll2007 

;|T 



.._ .,,i[P 

Blake et ail2007l1fegmark et al.l2006tlPercival et al.l2007h . The 



realization that BAO can be used as a standard ruler to map the 
low-r e dshift expansion history of the Universe ( Eisenstein & Hu 



1998 



2003 



Blake & Glazebrookl 120031; iHu & Haiman .2003; .Lindei 
Seo & EisensteinI 20031) . and thus to shed light on the 
nature of the mysterious dark energy (DE), has lead 



sev- 



eral groups to propose next generation redshift surveys that 
should cover even larger sky areas and go significantly deeper 
in redshift. A few examples of the proposed projects include 



' http://www.mso.anu.edu.au/2dFGRS/ 
^ http://www.sdss.org/ 
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BOSSB WFMO^i ADEPlH Hetde?(Q, and also already ongo- 
ing WiggleZQ project performed on Anglo- Australian telescope. 

In addition to these galaxy redshift surveys a fast de- 
velopment in the microwave detector technology has lead 
to the possibility of detecting tens of thousands of galaxy 
clusters over several thousand squa re degrees through ther- 
mal S unyaev-Zeldovich (SZ) effect ( Sunva ev & ZeldovichI 19721 
11980 ). Currently the most ambitious amongst these projects is 
the ongoing SP10 survey. Also there are plans in the X-ray com- 
munity to carry out next generation X-ray cluster survey with 
the yield of ~ 100, 000 galaxy clusters (eROSITA0 survey). With 
a further possibility of obtaining spectroscopic redshifts for the 
SZ/X-ray selected galaxy clusters through the optical follow-up, 
one can turn these large cluster samples to an extremely power- 
ful cosmological probe. 

In this paper we study the possible synergy between these 
galaxy and galaxy cluster surveys, particularly focusing on the 
two-point clustering measures. As the simplest approach one 
would analyze these surveys independently, i.e. concentrating 
separately on galaxy-galaxy and cluster-cluster pairs. However, 
this approach would neglect extra information provided by the 
cluster-galaxy cross-pairs. In fact, the study of th e cluster-galaxy 
cross-correlations has rather long history, e.g . PeebleslTl974l; 



ISeldner & Peebieslll977albt iLilie & Efstafliioij|[T988l Although 
rather sparse, cluster samples with strongly increased clustering 
strength as compared to the galaxies, had proven to be powerful 
tracers of the large-scale structure of the Universe. Initially the 
main interest in determining the cluster-galaxy cross-correlation 
function was driven by the desire to reduce the shot noise level 
inherent to these relatively sparse cluster samples. 

As the main task of this work we are going to determine the 
information gain provided by the cluster-galaxy cross-spectrum 
(Pc-g) over the information encoded in the galaxy-galaxy and 
cluster-cluster spectra (Pg-g and Pc-c, respectivel y). We carry 
out our modeling using a Halo Model (HM) (e.g. ISeliakI 12000. 
1200 U ICoorav & ShethI 12002 ) approach. To our knowledge this 
is the first time when HM has been used to model Pc-g- As such, 
we try to give a rather comprehensive explanation of the results 
obtained. We apply the derived formalism to a BOSS-like galaxy 
survey covering 10,000 deg^ and reaching redshifts of 0.7. For 
the cluster survey we assume that it covers the same volume 
as the galaxy survey and is simply mass-limited with the lower 
mass cut-ofFmJ.™ = 1.7 x IO'^^/i^'Mq, which is close to the mass 
limit obtainable by the SPT. In reality BOSS and SPT surveys 
do not overlap, but eROSITA, which has roughly the same effec- 
tive sensitivity (although the threshold mass in this case depends 
more strongly on redshift) is a full sky cluster survey, and thus 
has a complete overlap with BOSS. 

In our study we focus on two-point clustering statistics. Thus 
our approach is complementary to the number count analysis. In 
case of the galaxy cluster samples the complementarity of the 
number count and clusteri ng analysis has been investigated by 
iMajumdar & Mohj (120041) . 

Our paper is organized as follows: In Section 2 we present 
the formalism, which is applied to hypothetical galaxy and clus- 
ter cluster surveys in Section 3. There the main focus is to find 



out the gain in information once cluster-galaxy cross-pairs are 
included. Finally, Section 4 brings our conclusions. 

Throughout the paper we assume a fiducial cosmology to be 
a flat "concordance" ACDM model (IBahcall et al.lll999l) wifli 
Q„, = 0.27, = 0.045, h = 0.7, and erg = 0.85. All these 
parameter values, except somewhat larger value for erg, are in 
good agreement with the recent WMAP5 results jPunklev et al.l 
.2008.) . 

2. Halo Model power spectra. Covariances and 
Fisher matrices 

2.1. Power spectra in real space 

In this Subsection we present the results for the matter (m-m), 
galaxy (g-g), cluster (c-c), and cluster-galaxy (c-g) power spectra 
in real space within the framework of the Halo Model (HM). The 
corresponding results for the redshift-space are given in the fol- 
lowing Subsection. For a comprehensive review on HM see e.g. 
ICoorav & ShethI (I2002h . To the readers familiar with the HM we 
suggest to directly jump to Eqs. Q-®, which present a com- 
monly used formulation for the g-g power spectrum in real 
space. Our formalism for the c-c and c-g spectra in real space 
is new and given by Eqs. (fT3Tl-(fT8]l. Eqs. (|23])-(l34li provide anal- 
ogous formulae for the anisotropic spectra in redshift-space. 

In HM the power spectra can be expressed as a sum of the 
one-halo (Ih) and two-halo (2h) terms: 



P,- = Pi" + P 



2h 



(1) 



which provide the contributions of pairs of points inhabiting the 
same or separate dark matter halos, respectively. For our cases 
of interest the subscript x can be either m-m, g-g, c-c, or c-g. 

Following Seljak (2000) and lCoorav & ShethI (|2002|) the Ih 
and 2h terms for the matter power spectrum at redshift z can be 
expressed as: 



„(%) = J dmn(m\z)^^ ul,(k\m,z), 



Pftjklz) = bl(k\z)gHz)P'^(k\z = 0), 



bm(k\z) 



dm n(m\z)b(m\z) — 
P 



Um(k\m,z) ■ 



(2) 
(3) 
(4) 



http://cosmology.lbl.gov/BOSS/ 
http://arxiv.org/ftp/astro-ph/papers/05 10/05 10272.pdf 



In the above expressions n(m\z) is the mass function and b(m\z) 
halo bias par ameter at redshift z. We calcul ate these us i ng pre - 
scriptions by ISheth & TormenI ([l999) and ISheth et all (1200 lb . 
g(z) represents linear growth factor at redshift z normalized such 
that g{z = 0) = 1, and p = Q.,„Pc - f dm mn(m\z)b(m\z) 
gives the mean comoving matter density. u„i(k\m, z) is the nor- 
malized (i.e. lim Um(k\m, z) = 1) Fourier transform of the dark 

k — >i) 

matter density profile within a halo of m ass m at redshift z , 
which we take to be given by the the NEW jNavarro et al.ll 19971) 
form. The concentration parameter of the density profiles and its 
evolutio n with redshift is ass umed to take an analytic form as 
given in Bullock et a n (I200II) . To calculate linear power spec- 
trum P^"\k\z - 0) we assume power-law initial spectrum with 



spectral index ris - 1 and use fitting formulae for transfer func- 

tions T(k) as given bv lEisenstein & Hul (Il998l) . 

e.g. http://www.science.doe.gov/hep/hepap/feb2007/ hepap-bennett_feb0 7.ppfi e corresponding results for the g-g case can be expressed 
http://www.as.utexas.edu/hetdex/ as (see lCoorav & ShethI ( |2002|) for details): 

http://wigglez.swin.edu.au/ 

h,tp://pole.uchrcago.edu/ Cgikk) - f dm «(m|z) "/^'"^ |.,(^|m, z)P , 

http://www.mpe.mpg.de/erosita/MDD-6.pdi * * J n'^(z) 



(5) 
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Fig. 1. Power spectra and bias parameters in real- (left-hand panels) and redshift-space (right-hand panels). The short-dashed, 
dashed, and solid lines correspond to cluster, galaxy, and cross-spectrum, respectively. Where relevant we have shown contributions 
from the l/i and 2h terms separately. In lower panels a superscript "tot" is used to denote a total bias parameter, which includes both 
Ih and 2h contributions, in contrast to the bias parameters bg, be, and beg - -sjbgbc, which incorporate 2h contributions only. See 
the main text for the assumed values of the HOD and other parameters. 



P;tg(k\z) = bl(k\z)g\z)P''"(k\z = 0), 



r {Ng\m) 
bgik\z) = I dmn(m\z)b(m\z) _ ^ ^ Ug(k\m,z), 



ngiz) 



rigiz) = J dmn(m\z){Ng\m} . 



(6) 
(7) 

(8) 



These results follow directly from the corresponding results for 
the m-m case if one imagines the mass density field as composed 
of discrete particles, and thus can write m/p = N/n where is 
the number of particles in a halo with mass m and n is the av- 
erage number density of particles. To count the number of pairs 
correctly in this discrete case we have a second factorial moment 
(instead of simply a second moment) of the Halo Occupation 
Distribution (HOD) in Eq. (|5]l. In the above relations hg{z) is the 
mean number density of galaxies and Ug{k\m, z) is the normalized 
Fourier transform of the galaxy number density profile within a 
halo of mass m at redshift z. For Ug{k\m,z) we assume that it 
follows the corresponding density distribution for the dark mat- 
ter i.e. we take Ug{k\m,z) - u„,{k\m,z)- We further assume that 
the 1st and 2nd factorial moments of the HOD are given by the 



following parametric forms (e.g. ICooravll2004l) : 



{Ng\m) 







if m > mj?" 
if m < mj,™ , 



(9) 
(10) 

' otherwise . 



<A^^,(A^^ - l)|m) = l5\m){Ng\mf , 



This HOD model is consistent of assuming that for the high mass 
halos (m > 10'^ /i 'Mq) the HOD follows Poisson distribution 
whereas for the smaller masses the fluctuations in A^^ are sub- 
Poissonian. Thus in our case we have three free parameters de- 
scribing the HOD: a, niQ, and mj,°*. 

Here we have made an assumption that if there are any galax- 
ies in a halo then one of these is located in the central position. 
Thus for the \h term we have to deal with two types of pairs: 
(i) ones that involve the central object, and so probe directly the 
density profile of the halo (i.e. the power law index p - I'm Eq. 
dU), (ii) all the other pairs, which describe the convolution of 
the density profile with itself (i.e. p - 2). The detailed treatment 
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(see e.g. ICoorav & She"5ill2002h of the central objects leads to 
the following choice for p in Eq. ([5]): 



P = 



1 if {NgiNg - l)|m> < 1 

2 if {NgiNg - l)\,7t} > I 



(12) 



Following the derivations for the matte r and galaxy power 
spectra as given m ICoorav & ShethI (l2002h we obtain the cor- 
responding results for the c-c and c-g cases. The cluster power 
spectrum can be given as; 



Pl-Ak\z) = 

bc(k\z) = —— I dmn 
ndz) J 

ficiz) = J dmn(m\z) 



Pl'l,(k\z) = bl(k\z)gW'(k\z^O), 



(m\z)b(m\z) , 



(13) 
(14) 

(15) 



(16) 



and the cluster-galaxy cross-spectrum as: 



P'clgm = _ . . . r dmn{m\z){Ng\m)ug{k\m,z), (17) 
f^c\Z)ng\Z) J 

(18) 



Pf_g(k\z) = bMz)bg(k\z)gHz)P''"(k\z^O). 

It is clear that there is no Ih term for the cluster power spec- 
trum, [3 whereas the 2h term can be obtained from the corre- 
sponding g-g 2h term by replacing {Ng\m} with 1 (i.e. each dark 
matter halo hosts one cluster), hg(z) with hdz), and taking the 
normalized Fourier transform of the density profile (5-function 
in our case) equal to 1. [3 For calculating cluster power spectra 
we assume for simplicity that our cluster sample is mass-selected 
above some mass limit mj,™ which can be dependent on redshift. 

On small scales (i.e. relevant for the Ih term) cluster-galaxy 
cross correlations probe density profile around halo center i.e. 
in comparison to the g-g case, where power law index p in 
Eq. (|5]i can take value 2, here we have only a linear dependence 
on Ug{k\m,z). The number of pairs contributing to the \h term 
is simply given by the number of galaxies within each halo i.e. 
{Ng\m). 

In the above formulae we have used the notations b,n{k\z), 
bg(k\z), and bc{k\z) to emphasize the fact that these functions 
give the effective bias parameters for large scales (i.e. neglect- 
ing the contributions from l/i terms). From Eq. (fTSl l we see that 
the effective large-scale bias parameter bc-g{k\z) for the cross- 
spectrum is given by ^]bc{k\z)bg{k\z), as expected. 

Examples of the real-space spectra and bias parameters are 
given on the left-hand panels of Fig. [T] The results presented 
there correspond to the redshift z - 0.5. To calculate the cluster 
spectrum we have assumed a mass-limited sample of clusters 
above a threshold mass mj°* = 1. 7 xIO'^^/i^'Mq. The parameters 
for the galaxy HOD were fixed as follows: a - 0.9, m'™ - 
5 X 10'^ /j^'Mq, and mo - 4.2 x 10'^ /i^'Mq. For the motivation 



Here we have chosen to subtract the shot noise term, which effec- 
tively corresponds to the cluster number counts. Thus our clustering 
analysis is independent and complementary to the number count stud- 
ies. 

" To be more precise, even in Eq. (O one has to take Ug(k\m,z) = 1 
if (Ng\m) < 1, since as we assumed, if there is one galaxy per halo, it 
must be located at the central position. In reality, however, these details 
make a negligible difference. 



of these values see Sec. |3] In Fig. [T]the short-dashed, dashed, 
and solid lines correspond to cluster, galaxy, and cross-spectrum, 
respectively. Where relevant we have also shown contributions 
from the l/i and 2h terms separately. Here it is worth noticing a 
relatively bigger importance of the \h term in case of the cross- 
spectrum as compared to the galaxy-galaxy power spectrum. In 
the lower panel we have used a superscript "tot" to denote a total 
bias parameter, which includes both l/i and 2h contributions, in 
contrast to the bias parameters bg, be, (see Eqs. (fTt.lfTSl)) and 
beg — ^jbgbc, which incorporate 2h contributions only. For small 
values of the wavenumber k the total bias parameters deviate 
slightly from the results obtained by using 2h terms only due to 
the large-scale "leakage" of the \h term Pi 



2.2. Redshift-space power spectra 

In this Subsection we are going to generalize the previous re- 
sults to include the redshift-space distortions. Since in addition 
to the above-mentioned future redshift surveys there are several 
imaghm surveys planned (e.g. DE^l Pan-STARR^3 DUN^l 
LSSlQ, in Appendix|A]we also give a simple treatment in case 
only photometric redshifts a re available. H ere w e partially fol - 
low the formalism as given in lWhitd (1200 ll) and in lSeliakl(l200lb . 
For the results presented in this Subsection we assume a flat sky 
approximation. 

Small-scale redshift-space distortions ("Fingers-of-God" 
(FOG)) in coordinate space can be modelled as a convolution 
along the line of sight with a ID Gaussian kernel. As Gaussian 
stays Gaussian also in Fourier space the matter density contrast 
gets modified as: 



^ ■ nkcr„n) , 
where the damping factor ;F is given by: 

T(x) = exp -- 



(19) 



(20) 



Here fi is the cosine of the angle between the fine of sight di- 
rection f and wavevector k, i.e. yu = f ■ k. As in ICooravl (12004 . 
we assume that the one dimensional velocity dispersion of mat- 
ter inside a halo with mass m at redshift z, o-m(m,z), follows the 
scaling of the isothermal sphere model: 



o-m(m, z) = r ■ 



Gm 



2rvir(m, z) ' 



(21) 



where rvir is the virial radius of the halo and y is a free parameter, 
which in the following we keep fixed to 1, i.e. assume precisely 
the value for the isothermal sphere model. The corresponding 



It is clear that these additional Poissonian fluctuations on the largest 
scales are incompatible with the large-scale homogeneity requirement. 
This is a shortcoming of the HM that hasn't been solved satisfactorily so 
far. A possible way around it would be the use of compensated density 
profiles for the halos, i.e. no DC component in their Fourier transforms. 
Even then, in 2h terms one still has to use uncompensated profiles, since 
convolving with the profile with zero DC component would damp out 
the large-scale power completely. However, for the scales of interest in 
this paper the deviations of the bias parameters are rather small, and 
thus we do not make an attempt to correct for this shortcoming. 

https://www.darkenergysurvey.org/ 

http://pan-starrs.ifa.hawaii.edu/public/ 

http://www.dune-mission.net/ 

http://www.lsst.org/ 
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velocity dispersion of galaxies is taken to follow the same rela- 
tion. 

There are additional redshift-space distortions on large 
scales due to the coherent inflows of matter towards massive 
accretion centers. This leads to the m odification of the density 
fluctuations in Fourier space such that (lKaisedll987h : 



(22) 



where / represents the logarithmic derivative of the growth fac- 
tor g{z), i.e. / = d log ^/d logo, and superscript x can be either 
m (for matter), c (clusters), or g (galaxies). 

Taking into account both the small and large-scale distor- 
tions we can write down the results for the anisotropic power 
spectra as follows; 



n(k,n\z) 



(23) 



dm n{m\z) 



{NgiNg - \)\m) 



nliz) 



Ug(k\m,z)f 



rP[kcT,„(m,z)fi] 
1 



nc(z)ng(z) 



= bcik,fj.\z) + f(z)b„(k,fi\z)ir 
= bg(k,fi\z) + f(z)b„(k,fi\z)i/ 



■ T[ko-Jtn,z)ij] , 
PlUk,n\z) = bl(k,ti\zfg\z)P''^(k\z = 0), 
Pf_g{k,n\z) = fo^(fe,//k)Vfe)^""(% = 0), 

pf_g{k,n\z) = fo^;^(^,//k)Vfe)^""(^k = o), 
K{k,^l\z) 

b\{k,n\z) 
Kg(k,lAz) 

bc(k,M\z) 
bg{k,n\z) 

b,„{k,ii\z) 



^blik,fi\z)bl(k,^i\z) 



ndz) 



dm n(m\z)b(m\z) .. 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 

(31) 
(32) 



m!"*(z) 



dm n(m\z)b(m\z) 



{Ng\m} 
ngiz) 



Ug(k\m,z) 



T[ka-,n{m,z)n] , 



dm n(m\z)b(m\z) — u,„{k\m, z) 



P 



T [kcr„,{m, z)iu] 



(34) 



Examples of the spectroscopic redshift-space spectra and 
bias parameters are given on the right-hand panels of Fig. [1] 
The HOD and other parameters are exactly the same as on the 
left-hand panels (see the end of Subsec. IZTI for details), which 
show the corresponding results for the real-space. As before, a 
superscript "tot" on lower panel denotes a total bias parameter, 
which includes both Ih and 2h contributions, in contrast to the 
bias parameters b^, b^, and b'^g (see Eqs. (I29b.(l30]l.(l3ni). which 
incorporate 2h contributions only. The spectra and bias param- 
eters presented here are the angle-averaged versions. Note the 
expected large-scale boost and small-scale damping of power as 
compared to the equivalent real-space results. 

2.3. Cosmological distortion 

The cosmological distortion (lAlcock & Paczvn ski"l979^ arises 
due to the simple fact that conversion of the observed redshifts to 



comoving distances requires the specification of the cosmologi- 
cal model. If this cosmology differs from the true one, we are left 
with additional distortion of distances along and perpendicular to 
the line of sight. In general, the spatial power spectrum measure- 
ments, in contrast to the angular spectra, are model dependent 
i.e. along with the measurements of the 3D power spectrum one 
always has to specify the so-called fiducial model used to ana- 
lyze the data. In principle, for each of the assumed cosmological 
models one should redo the full power spectrum analysis to ac- 
commodate different distance-redshift relation. However, there 
is an easier way around: one can find an approximate analytical 
transformation that describes how the model spectrum should 
look like under the distance-redshift relation given by the fidu- 
cial model. This transformation for the power spectrum is dif- 
ferent along and perpendicular to the line of sight, and also is 
dependent on redshift. For more discussion on these issues see 
e.g. lHutsil(l2006bh . The observed spectrum P ^ under the assumed 
distance-redshift relation of the fiducial c osmological m odel is 
related to the true spectrum as follows (lHutsill2006bh : 



(24) P,(fc||,fc^|z): 



1 



C||(z) ■ cl(z) 



k^ 



c\\(z) Ci(z) 



(35) 



where the distortion parameters along and perpendicular to the 
line of sight are given as: 



c\\(z) 
cAz) 



H^%z) 
H(z) 
dAz) 
df{z) ■ 



(36) 
(37) 



Here H{z) denotes the Hubble parameter and t/j.(z) is the co- 
moving angular diameter distance. Superscript '^'^ refers to the 
fiducial model. Here and in the following we use a tilde on top 
of Pf to denote theoretical spectrum "transformed to the refer- 
ence frame of the fiducial cosmology". As we use the spectra that 
have the dimensions of volume an extra division by c\\{z) ■ c\{z) 
occurs due to the transformation of the volume elements. The 
same transformation in terms of k - (A;^ -i- A;^)'^^ and n - k^/k 
can be expressed as: 



(33) PAk,f^\z) = 



1 



C||(z) ■ cl(z) 



■p. 



ai^i, z) 



k. 



1 



c±(z) x(z)a(fi,z) 



where 



x(z) = 



c\\(z) 
cAz) 



a(fi,z) = Jl 



x\z) 



(38) 

(39) 
(40) 



2.4. Covariances and Fislier matrices 

Having a HM description for the various spectra we can go on to 
calculate covariances and Fisher matrices. We give all the nec- 
essary details in AppendixiBl 



3. Application of the formalism 

In this Section we are going to apply the machinery presented in 
the previous section to the hypothetical ideal galaxy and galaxy 
cluster redshift surveys. For the galaxy survey we assume spec- 
ifications similar to the planned BOSS survey: 10, 000 deg^ sky 
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Fig. 2. Comoving number density (upper panel) and the prod- 
uct of the real-space bias parameter and linear growth factor 
b{z)giz) (lower panel) for galaxies (solid lines) and galaxy clus- 
ters (dashed lines) as a function of redshift. 



coverage, redshift range z = 0. 1 - 0.7, and the comoving num- 
ber density of galaxies ~ 3 x 10""* Z;-' Mpc""'. Thus the volume 
of the survey is ~ 5.7 h'^ Gpc^ with the total number of galax- 
ies reaching ~ 1,700,000. As mentioned in Sec.|2]we assume 
a simple HOD given by Eq. (|9|l, i.e. we have three free parame- 
ters: a, mo, and ;7i'°*. For the parameter a we choose a fiducial 
value of 0.9, which is typical for the red galaxie s according to the 
semi- analytic galaxy formation models (see e.g. lCoorav & ShethI 
I2OO2IFI . Looking at the equations given in the previous section 
it is evident that power spectra do not depend on niQ. Thus af- 
ter fixing the value for a the fiducial value for the parameter 



m;°* can be chosen such as to give a clustering strength typical 
of the luminous red galaxies targeted by BOSS, i.e. the product 
of the bias and linear growth factor b{z)giz) ~ 2. This gives us 
^low ~ 5 X 10'^ Iz'^Mq. Once a and m^°" are fixed the fiducial 
value for niQ is obtained through the number density constraint, 
giving mo ^ 4.2 x lO'^/i^'M©. The total survey volume is di- 
vided into six redshift bins with a width Az = 0.1. 



" However, from the observational data se veral authors have obtained 
values for a which are somewhat higher: ICollister & Lahavl ( l2005h 
found a = 1.05 for the 2dFGRS red galaxies jBlake et alj ilOtM found 
a for the LRGs to be evolving with redshift, from 1.57 to 1.80. 



The galaxy cluster survey is assumed to cover exactly the 
same volume. For simplicity we model our cluster survey as a 
mass limited survey above a threshold mass of 1.7 x 10'"* /i^'M©. 
This is very similar to the mass sensitivity obtainable by the SPT 
survey. Of course one could model the cluster selection (e.g. SZ 
or X-ray selection) more realistically taking the threshold mass 
to be an appropriate function of redshift. However, for clarity 
we choose to keep things as simple as possible. As an addi- 
tional prior we assume that the threshold mass can be calibrated 
with an accuracy of 10%. It turns out that the information in 
the cluster-galaxy cross-spectrum itself helps in determining the 
threshold mass with an accuracy of ~ 4%, and thus our results 
are rather insensitive to this assumed prior 

In addition to the HM parameters we have to fix a fiducial 
cosmological model. This we take to be a flat ACDM model 
with 0,„ = 0.27, Q,, = 0.045, h = 0.7, and os = 0.85. We 
will focus only on spatially flat models in our Fisher matrix 
calculations. For the dark energy equation of state we use a 
common wq - Wa parametrization: w(a) - wo + (1 - a)wa 
(a = 1/(1 + z)). So our parameter vector has nine compo- 
nents in total: = {Q.„„Q.i„ h,o-s,wo,Wa,\nm^°'^ ,a,lnm^c^) 
with the following fiducial values: 0*''' = 
(0.27,0.045,0.7,0.85, -l.,0.,ln(5 X 10i-),0.9,ln(1.7 x lO'^)). 

Having fixed the parameters to those fiducial values we end 
up with a total of ~ 1,700,000 galaxies and ~ 25,000 galaxy 
clusters. The redshift dependence of the comoving number den- 
sity and the product of the real-space linear bias parameter b{z) 
and the linear growth factor g{z) (i.e. the linear bias parameter 
with respect to the z - matter power spectrum) is shown on 
the upper and lower panels of Fig. |2] respectively. In redshift- 
space the product b{z)g(z) gets boosted by a factor of 1.12 and 
reaches ~ 2.1 at redshift z = 0.1 as compared to the correspond- 
ing real-space value of 1 .85. 

The important ingredients of the Fisher matrix calculations 
are the weight factors (w*, w^, w^"*) and the logarithmic deriva- 
tives of the power spectra (see Eqs. ( |B.5t ,( |B.6l l,( IB.7t ,( |B?8] l). 
The angle-averaged weight factors as functions of the comov- 
ing wavenumber k and redshift z for the galaxies, clusters, and 
cluster-galaxy cross-pairs are shown as upper, middle, and lower 
panel of Fig. |3] respectively. There the discrete contour lines 
start from the value of 0.5 and increase with a constant step 
of 0.05. As expected, we see that the number density of galax- 
ies (~ 3 X 10"^ /j^ Mpc"^) is sufficiently high to ensure a well- 
sampled density field (i.e. weight factors close to 1) for all of the 
wavenumbers and redshifts of interest. This is not the case for 
the galaxy clusters for which the high shot-noise contribution 
reduces the weights significantly, reaching the value of 0.5 at 
large scales only for the lowest redshifts. For the cluster-galaxy 
cross-pairs the situation is significantly better: here the weight 
factor is above the value of 0.5 for the most of the k - z plane. 

The logarithmic derivatives of the galaxy, cluster, and cross- 
spectra are given in Fig.|4]as dashed, dotted, and solid lines, re- 
spectively. We see that for most of the parameters the derivatives 
are rather similar, which leads to almost identical degenerate pa- 
rameter combinations. However, there are significant differences 
in logarithmic derivatives involving and HOD parameters 
a and lnm|,°*, and thus for these one would expect tighter con- 
straints once the extra information from the cross-pairs is in- 



Note the difference in sign of the logarithmic derivatives of the 
galaxy and cluster power spectra in this case. By increasing erg the am- 
plitude of the cluster power spectrum actually drops. This owes to the 
fact that the higher value of erg leads to the rapid increase of the number 
of low mass clusters which are less weakly biased. 
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eluded. This is indeed the case as can be seen from Fig. |5] and 
Table [T] In Fig. |5]we have plotted the 2D error ellipses (1-cr) 
and also ID Gaussian distributions for the eight free parameters 
involved. There the largest ellipses correspond to the constraints 
from the cross-pairs, middle ones to the galaxy-galaxy pairs, and 
the smallest to the combined (g- g) + {c- g) case. In Table[T]we 
give the figures-of-merit (FOMs), defined as the inverse areas of 
the 1-cr error ellipses, for several free parameter combinations. 
The diagonal entries (e.g. Q„,-Q„,) in this table give the inverse 
of the corresponding parameter's 1-cr error interval, instead. In 
each cell the numbers from top to bottom correspond to c - g, 
8-8^ {8-8) + {c-g), and the ratio of (g-g) + (c-g) snAg-g 
FOMs. In Table [1] we have decided not to show an additional 
gain once c - c information is included, as due to the sparse- 
ness of the cluster sample this is relatively modest compared to 
the accuracy boost brought by the inclusion of the cluster-galaxy 
cross-pairs. As already mentioned, the biggest boost in accuracy 
once cross-pairs are included occurs for the parameter combina- 
tions involving erg, a, and Inrnj,"* as at least one of the param- 
eters, e.g. the error ellipse for the ct^-wq case gets reduced by a 
factor of ~ 2.4, etc. Since the g- g spectrum does not depend on 
the ninth parameter In m|°* we have decided not to include it in 
Fig. |5] and Table[T] We only note that even without the assumed 
10% prior for mj,°* the information encoded in the c - g pairs is 
able to self-calibrate mj.°™ to an accuracy of ~ 4%, and thus our 
results in Fig.|5]and Table [1] change only slightly once this extra 
prior is not applied. 

In Fig. |6] we show in greater detail the improvement on the 
m'°*-Q; plane brought by the inclusion of the information from 
the cross-pairs. Although the constraints in the c- g case are rel- 
atively weak a different degeneracy direction as compared to the 
g—g case still helps in shrinking ^g-g error ellipse by a factor of 



~ 2 (see Table[T]). Also, as one of the main goals of these future 
redshift surveys is to put constraints on the properties of dark 
energy. Fig. |7] presents the performance on the vvq-Wo plane in 
greater detail. The dark energy FOMs we obtain are 26.9, 38.5, 
28.8, and 40.4 for the g-g,{g-g) + {c- g), (2- g) + (c- c), and 
(g - g) + (c - c) + (c - g) cases, respectively^ Here we see that 
the extra information from the cross-pairs leads to a relatively 
modest gain: error ellipse is compressed by a factor of ~ 1 .4 as 
compared to the g- g case taken separately. We also see that the 
additional gain provided by the inclusion of the cluster-cluster 
pairs is significantly smaller than that. This is also the case for 
the rest of the parameters, which is why we have not included 
c - c case in Fig. |5] and Table [T] 

In our Fisher matrix calculations we have focused only on 
a wavenumber interval k = 0.01 - 0.2/iMpc"'. Due to forma- 
tion of the nonlinear structures the Gaussianity assumptions is 
not adequate for wavenumbers larger than ~ 0.2/!Mpc"'. At 
those small scales the information is transferred out of the two- 
point function to the higher order moments, increasing the vari- 
ance of the power spectrum above it's Gaussian expectation, 
i.e. to calculate the variance one has to in clude also a contri- 
bution from the trispectrum (see e.g. IScoc cimarro et al. 199^. 
Moreove r, there ha ve been a few studies (e.g. lRimes & Hamilton! 
2006; Ne vrinck & S zapudi 2007) that indicate that there is not 
much information left in the power spectrum at those scales. 
As we have chosen to divide our survey volume into rela- 
tively narrow redshift shells it is clear that we are not able to 
probe the very large Fourier modes. However, the modes with 

" For the flat models with a scalar spectral index = 1 (as was 
assumed so far) by adding extra information obtainable from the CMB 
measurements by the Planck satellite these FOMs increase to 81.9, 125, 
86.1, and 135, respectively. 
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Fig. 5. 2D error ellipses (1-cr) and ID Gaussian distributions for eight of the free model parameters. The largest ellipses correspond 
to the constraints from the cluster-galaxy cross-pairs, middle ones to the galaxy-galaxy pairs, and the smallest to the combined 
- ^) + (c - case. 



wavenumbers above 0.01 /zMpc" are still relatively well sam- 
pled. We have checked that our results are rather insensitive to 
lowering the wavenumber boundary below it's assumed value of 
k - 0.01 h Mpc"' . For the c - g case this can be seen by looking 
at Fig. [8] where we have shown how Fisher information (for the 
diagonal matrix entries only) is distributed on the wavenumber- 
redshift plane, i.e. the total Fisher information is given by the 
integral over this plane. There the discrete contour lines corre- 
spond to the isocontours enclosing 50, 70, and 90% of the total 
Fisher information. One can see that for most of the parameters 
the surface drops off rather rapidly towards large scales and thus 
90% contour does not reach k - 0.01 h Mpc ' line. For the ^ 
and c - c case the picture is qualitatively similar 



approach is complementary to the number count analysis, which 
in its usual form assumes validity of the Poissonian statistics. PI 

To model the cluster-galaxy cross-spectrum we have used 
the Halo Model framework. We have carried out a Fisher matrix 
analysis for a BOSS-like galaxy redshift survey targeting lumi- 
nous red galaxies over 10,000 deg2 of sky and over the redshift 
range of z = 0.1 - 0.7, and a hypothetical mass-limited cluster 
redshift survey (with a threshold mass m[°* = 1 .7 x lO'"* /i^'Mq) 
over the same volume. For simplicity we have assumed spatially 
flat models with the usual wq - Wa parametrization for the dark 
energy equation of state. For the HOD we have used a simple de- 
scription given by Eqs. (|9]l,(fT0ll,(fTni. Since clustering does not 
depend on mo we are left with two HOD parameters: a and 
7M^°*. Adding to these m[°* and six cosmological parameters: 
Qm,^^/),/',cr8,wo, and Wa, we have in total nine free parameters. 



4. Conclusions and Discussion 

There are several wide field galaxy and cluster surveys planned 
for the nearest future. In the simplest approach one would ana- 
lyze these independently, thus neglecting the extra information 
provided by the cluster-galaxy cross-pairs. In this paper we have 
focused on the possible synergy between these surveys by inves- 
tigating the amount of information encoded in the cross-pairs. 
Since in our study we focus on two-point clustering statistics our 



In case of the galaxy cluster samples the complementarity of 
the number count and cl ustering analysis has been investigated by 
iMaiumdar & Moh3 ( l2004h . For example, these authors illustrate that 
combining number counts and c — c power spectrum helps in reduc- 
ing the error bar on erg by a factor of two for an SPT-like survey. This 
is similar to the improvement on erg we obtain by combining c - g with 
g-g- Due to the sparseness of the cluster sample the additional inclusion 
of the c - c power leads to only mildly stronger constraints. To evaluate 
the absolute error on a% one should combine our pair statistics with the 
cluster number counts, which is beyond the scope of the current work. 
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Table 1. Figures-of-merit (FOMs), defined as the inverse areas 
of the 1-cr error ellipses, for most of the free parameter combina- 
tions. The diagonal entries (e.g. 0„,-Q„,) give the inverse of the 
corresponding parameter's 1-cr error interval, instead. In each 
cell the numbers from top to bottom correspond to c - g, g - g, 
ig-g) + {c- g), and the ratio of (g-g) + (c- g) and g-g FOMs. 



0"8 



Wo 



InrnJ,™ 



165 
232 
290 
1.25 



n,, 



5610 
12,300 
19,100 
1.55 



245 
384 
477 
1.24 



673 
1510 
2570 
1.70 



1550 
4020 
6580 
1.64 



26.9 
43.7 
58.7 

1.34 



506 
1420 
3450 
2.43 



736 
2470 
5700 
2.31 



87.2 
328 
758 
2.31 



21.6 
44.2 
85.6 
1.94 



Wo 



244 
523 
797 
1.53 



328 
827 
1250 
1.51 



34.0 
88.7 
144 
1.62 



30.8 
92.8 
226 
2.44 



9.13 
14.6 
17.6 
1.21 



60.4 
129 
197 
1.53 



73.2 
189 
282 
1.50 



7.81 
20.9 
33.5 
1.60 



6.88 
22.2 
53.2 
2.39 



9.86 
26.9 
38.5 
1.43 



2.09 
3.45 
4.12 
1.19 



In m'': 



52.5 
546 
998 
1.83 



78.3 
903 
1600 
1.77 



9.49 
107 
199 
1.86 



16.9 
163 
344 
2.12 



3.00 
34.1 
59.9 
1.76 



0.675 
8.06 
13.9 
1.73 



2.27 
16.6 
24.2 
1.46 



one cannot reliably perform profile fitting on object by object 
basis, and thus probably has to rely on stacking methods, which 
also lead to some sort of a weighted profile. Thus these two ap- 
proaches should give essentially the same results. 

In our analysis we have assumed that the threshold cluster 
mass m[.°™ can be calibrated with an accuracy of 10%. We find 
that even without the assumed 10% prior for m[°* the informa- 
tion encoded in the c - g pairs is able to self-calibrate mj.°™ to an 



612 accuracy of ~ 4%, and thus our results are rather insensitive to 
I67(lhis additional prior On the other hand, treating m[°* as a known 
347(i|uantity and fixing it to the fiducial value of 1.7 x IO'^^/i^'Mq 
^•*^^e obtain ~ 5 times tighter constraint on cr^. Similarly, all the 
error ellipses involving erg as one of the parameters get reduced 



970 
2951 



5850^^ factors of 6 - 10 in this case, the most remarkable being a fac- 
j ^^or of ~ 10 reduction of the error ellipse on the crg-m'°* plane. 

127 For the other parameters the precise knowlwdge of ot[°* leads to 
372 much milder improvements, e.g. m^°'^-a and wo-Wa error ellipses 
754 now decrease by factors of ~ 3 and ~ 1.5, respectively. 
2.03 In our analysis we have simply added Fisher matrices cor- 



87. 
581 
126q 



Responding to the galaxy-galaxy and cluster-galaxy pairs. This 
' is a valid approach once the density field is not over-sampled, 
■^^f that is not the case, the extra pairs do not help in increasing 



33.6, 



■tbe information, i.e. then we cannot simply add Fisher matrices. 
However, in our case the Fisher information comes mostly from 
2n the larger wavenumbers as there are simply more high-fc modes 
2.0(Fvailable (Fig. |8] shows this for the c - g case). As can be seen 
7.7()from Fig.|3] for large wavenumbers our hypothetical surveys are 
25.0surely not over-sampled, and thus simply adding the Fisher ma- 
49.3trices seems to be a valid approximation. 

As a final comment we point out that the formalism as pre- 



^l-^sented in Sec. |2] can also be applied to the photometric red- 
320 shift surveys such as DBS, Pan-STARRS, DUNE, and LSST. In 



^52|his case including cluster-galaxy pairs might be very beneficial 
since the photo-z errors for the clusters are usually significantly 
'smaller than that for the typical galaxies. 



26.6 
52.1 
86.4 
1.66. 



As on small scales the cluster-galaxy cross-spectrum probes 
directly density profiles of the halos, instead of the profile con- 
volved with itself, as is the case for the galaxy power spectrum, 
we are sensitive to a different combination of the HOD param- 
eters. This leads to stronger constraints on the HOD parameters 
a and mj,™ as compared to the constraints obtainable from the 
galaxy data taken separately, e.g. the inclusion of the cross-pairs 
leads to a factor of ~ 2 compression of the error ellipses on the 
m'°*-a' plane. Also, the extra information in the cross-pairs leads 
to a factor of ~ 2 tighter constraint on erg . The best performance 
is obtained on crg-f2,„, crg-wo, and crg-Wa planes where the error 
ellipses get reduced by a factor of ~ 2.4. 

For the dark energy equation of state parameters wq and w,, 
the improvement is somewhat weaker: inclusion of the cluster- 
galaxy cross-pairs leads to a factor of ~ 1 .4 increase in the dark 
energy figure-of-merit as compared to the galaxy-galaxy case 
taken separately. For the rest of the results see Fig. |5] and Table 
□ 

It is worth mentioning that the small scale information in the 
one-halo term of the cluster-galaxy power spectrum is essentially 
equivalent to the information gained by directly fitting for the 
galaxy distribution within the halos. P'.'!^ is just a Fourier trans- 
form of the weighted mean of the halo density profiles. From 
the rather sparse data provided by the assumed redshift surveys 
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Appendix A: Power spectra in photometric 
redshift-space 

As there will be several large imaging surveys performed in the 
future (e.g. DES, Pan-STARRS, DUNE, LSST), we present here 
a minimal extension of the formalism for the case when only 
photometric redshifts are available. For simplicity we assume 
that the conditional probability distribution for the spectroscopic 
redshift given the photometric one P(Zspeckphoto) is a Gaussian 
with a mean value of Zphoto (i-S- we assume zero bias) and with 
dispersion Sz- This leads to the extra damping factors for the 
density contrast, which have the same form as in Eq. ( fT9] l with 
the ID spatial smoothing scale cr given as cr = We denote 
the corresponding smoothing scales for clusters and galaxies as 
crf^°^° and crf^°^°, respectively. As for any realistic case photo- 

z errors change with redshift, crf^°^° and crf^°^° are modeled as 
functions of z- 

Thus for the dark matter we have an effective 
smoothing scale cr„,(m,z), for the galaxies o-g(m,z) = 

■\lio-m{m,z))^ + (cr^^°'°(z))^ (FOG + photo-z errors), and 
for the clusters o-dz) = crf^°^°(z) (photo-z errors only). We 
also define the quantity o-c-g(m,z) = yjio-dz))^ + {a-g{m,zS9- - 



(crm(OT,z))2 H- (cr5!'""°(z))2 -H (o-^""'"(z))2. Having introduced 
these new smoothing scales we can obtain results valid for 
the photometric redshift-space if we substitute cr„,(m,z) with 
o-f,(m,z) in Eqs. ( l24l l and (|33] |. cr„j(m,z) with crc-g(m,z) in Eq. 
( |25l l, and include an additional factor of f [ka-c(z)lj] in Eq. ( [32] i. 



where V is the total volume of the survey. Here the second equal- 
ity for VeSf{k) is valid only if the selection function of the clus- 
ter/galaxy sample, n,i(r), is independent of position. 

Following the l ines of thought in Appendix B of 
iFeldman et al.l (1 19941) . which lead to the result given in Eq. 
(IB. Il l we find that the similar result for the cluster-galaxy cross- 
spectrum can be expressed as: 



sr - 



Pi-g{k„) _^ Pc-c(kn)Pg-g(kn) 



v„v 



v„yv(fc„) 



(B.3) 



where the effective volume V^^^ik): 



v:,^k) 



h,(r)P,^,(k) 



1 + nc(r)Pc- 
hcPc-c(k) 



c(k) 



ng{r)PgJk) 



1 + ng{r)Pg-g{k) 



dr 



1 + n,Pc-c{k) 



^gPg-gik) 



1 + figPgJk) 



y. 



(B.4) 



The first term in Eq. ( IB.3I ) is the analog of Eq. dB.ll ) with the 
difference of the missing factor of two and also that here in- 
stead of the effective volume we have a survey volume V, which 
is due to the fact that in case of the cross-spectrum we do not 
have a shot-noise term raising from self-pairs. To fin d the cor- 
respon ding Fisher matrices we follow the derivation in lTegmarfl 
(fT997l) .In our calculations we divide the survey volume into rel- 
atively narrow redshift bins. This allows us to approximate the 
spatial number densities of objects as constant quantities within 
the bins, i.e. we can use the second equalities of Eqs. ( |B.2| i and 
(El. 

The Fisher matrices for the anisotropic cluster/galaxy spectra 
can be given as: 



A "max 

'"F,/z„) - ^ ■ / ''^ / '^^^^'^k,n\z„) 
d\nPAk,n\z„) din P,(k,ii\z„) 



w''(k,fl\Zn) = 



50, 

hJz„)Px(k,fi\Zn) 



1 + n,(z„)Px(k,iJ\z„) 



d&j 

v^(Mk«) 

V(Zn) 



(B.5) 



(B.6) 



Here z„ denotes the central redshift and V(zn) the comoving vol- 
ume corresponding to the n-th redshift bin. 

The corresponding results for the cluster-galaxy cross- 
spectrum are given as: 



Appendix B: Covariances and Fisher matrices 

Under the assumption of Gaussianity the p ower spect rum covari- 
ance matrix can be expressed as dFeldman et al.lll994; ,Tegmari3 
119971) : 



, 2Pl(k„) ^ 



where Px(k„) denotes the power spectrum of clusters/galaxies 
for the wavenumber bin centered at wavenumber k„, V„ = 
Ank^Akn / (2jt)^ {Ak„ - width of the wavenumber bin) is the vol- 
ume of the shell in Fourier space, and the effective volume V-^g{k) 
is given as: 



'-'FijiZn) 



^ ■ jdf^j dkk^w'-Hk,fi\zn) 
8 In Pc-g{k,iJ.\z„) d In Pc-g(k,fi\z„) 



d&i 



d@j 



(B.l) w'-Hk,tl\Zn) = 



, [l+«,fe„)f.-.(<:,Ai|z„)][l+«g(z„)f„(<:,//k„)] 
«r(z„)ns(z„)/'?-,(<:,/i|z») 



(B.7) 
(B.8) 



The total Fisher matrices are then the sums over the ones corre- 
sponding to different redshift bins: 



'F^j^J]^F^j(z„). 



(B.9) 



v:Ak) 



nAr)Px(k) 



1 + hx(r)PAk) 



fixPxik) 



1 + hxPxik) 



y. 



(B.2) 



In the above relations 0, represents the /-th component of the 
parameter vector 0, which in our case consists of cosmological 
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plus HM parameters. According to the Cramer-Rao inequality 
there does not exist an unbiased method that can measure the 
i-th parameter 0, with error bars less than 1 / -\^. If the other 
parameters are not known but estimated from the data as well, 

the minimum achievable error is given by ^F^' instead. 




11 m i,» irf *i m- m 




m m m- if 




li m m m m m 



i: 

Fig. 3. The angle-averaged weight factors (upper panel), w'^ 
(middle panel), and w'^^^ (lower panel) as functions of the co- 
moving wavenumber k and redshift z. The discrete contour lines 
start from the value of 0.5 and increase with a constant step of 
0.05. 
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Fig. 6. 1 -cr error contours on the m^"™-^ plane. Solid, dashed, and 
short-dashed lines correspond to the g-g, c-g, and (g-g)+ic-g) 
cases, respectively. 
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Fig. 7. 1-<T error contours on the wo-Wa plane. Solid, dashed, 
short-dashed, and dotted lines correspond to the g - g, (g - g) + 
(c - g), (g- g) + ic- c), and (g - g) + (c - c) + (c - g) cases, 
respectively. 
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Fig. 8. Distribution of the c-g Fisher information (for the di- 
agonal matrix entries only) on the wavenumber-redshift plane, 
i.e. the total Fisher information is given by the integral over this 
plane.The discrete contour Unes correspond to the isocontours 
enclosing 50, 70, and 90% of the total Fisher information. 



